THERMINATOR: 
THERMal heavy-IoN generATORt 



Y^-, Adam Kisiel", Tomasz Taluc", Wojciech Broniowski and 

O ' Wojciech Florkowski "^'^ 

o ■ 

Faculty of Physics, Warsaw University of Technology, PL-00661 Warsaw, 

Poland 

C/^ ! ^ The H. Niewodniczanski Institute of Nuclear Physics, 

\ Polish Academy of Sciences, PL-31342 Cracow, Poland 

^ Institute of Physics, Swigtokrzyska Academy, PL- 25406 Kielce, Poland 

(N 
> 

o 

■ THERMINATOR is a Monte Carlo event generator designed for study- 

. ing of particle production in relativistic heavy-ion collisions performed 

at such experimental facilities as the SPS, RHIC, or LHC. The pro- 
gram implements thermal models of particle production with single 
freeze-out. It performs the following tasks: 1) generation of stable par- 
^ , tides and unstable resonances at the chosen freeze-out hypersurface 

^ ' with the local phase-space density of particles given by the statistical 

C ■ distribution factors, 2) subsequent space-time evolution and decays 

of hadronic resonances in cascades, 3) calculation of the transverse- 
\ momentum spectra and numerous other observables related to the 

j_j ■ space-time evolution. The geometry of the freeze-out hypersurface and 
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Abstract 



the collective velocity of expansion may be chosen from two success- 
ful models, the Cracow single- freeze-out model and the Blast-Wave 
model. All particles from the Particle Data Tables are used. The 
code is written in the object-oriented C++ language and complies to 
the standards of the ROOT environment. 
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PROGRAM SUMMARY 



Manuscript title: Therminator - THERMal heavy loN generATOR 

Authors: Adam Kisiel, Tomasz Taluc, Wojciech Broniowski, Wojciech 
Florkowski 

Program title: THERMINATOR, April 2005, version 1.0; 

RAM required to execute with typical data: 50 Mbytes 
Number of processors used: 1 

Computer(s) for which the program has been designed: 

PC, Pentium III,IV, or Athlon, 512 MB RAM not hardware dependent; 

(any computer with the C++ compiler and the ROOT environment [1]) 

Operating system(s) for which the program has been designed: 
Linux: Mandrake 9.0, Debian 3.0, SuSE 9.0, Red Hat FEDORA 3, etc., 
Windows XP with Cygwin ver. 1.5.13-1 and gcc ver. 3.3.3 (cygwin special); 
- NOT SYSTEM DEPENDENT 

External routines/libraries used: 
ROOT ver. 4.02.00; 

Programming language: C++ 

Size of the package: (324 KB directory 40 KB compressed distribution archive), 
without the ROOT libraries (see |http: / / root.cern.ch for details on the ROOT 
[1] requirements). The output files created by the code need 1.1 GB for each 
500 events. 

Distribution format: tar gzip file 

Number of lines in distributed program, including test data, etc.: 6516 

Keywords: Monte Carlo generator, relativistic heavy-ion collisions, thermal 
models, particle production, transverse-momentum spectra 

PACS: 25.75.-q, 25.75.Dw, 25.75.Gz, 25.75.Ld 

Nature of the physical problem: 

Statistical models have proved to be very useful in the description of soft 
physics in relativistic heavy-ion collisions [2]. In particular, with a few phys- 
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ical input parameters, such as the temperature, chemical potentials, and 
velocity of the collective flow, the models reproduce the observed particle 
abundances [3-15], the transverse-momentum spectra [16], balance func- 
tions [17,18], or the elliptic flow [19,20] in both non-strange and strange 
sectors. The key element of the approach is the inclusion of the complete list 
of hadronic resonances, which at the rather high temperature at freeze-out, 
~ 165 MeV, contribute very significantly to the observed quantities. Their 
two- and three-body decays, taken from the tables, proceed in cascades, ulti- 
mately producing the stable particles observed in detectors. At the moment 
there exist several codes to compute the abundances of particles (the publicly 
available programs for this purpose are SHARE [21] and THERMUS [22]), which 
is a rather simple task, since the abundances are insensitive to the geometry 
of the fireball and its expansion. On the other hand, the calculation of the 
transverse-momentum spectra of particles is much more complicated due to 
the sensitivity to these phenomena. THERMINATOR deals with this problem, 
offering the full information on the space-time positions and momenta of the 
produced particles. As a result, the program allows to compute very effi- 
ciently the transverse-momentum spectra of identified particles and examine 
implications of the assumed expansion model. THERMINATOR allows easily 
for the departure from symmetries typically assumed in other approaches. 
This opens the possibihty to study the dependence of physical quantities 
on rapidity and the azimuthal angle. The contribution of the resonances to 
various observables may be traced conveniently, and their role in the statis- 
tical approach may be verified. As a Monte Carlo event generator written in 
the object-oriented C++ language in the ROOT [1] environment, THERMINATOR 
can be straightforwardly interfaced to the standard software routinely used 
in the data analysis for relativistic heavy-ion colliders, such as SPS, RHIC, 
and, in the future, LHC. In this way the inclusion of experimental acceptance, 
kinematic cuts, or interfacing with other programs poses no difficulty. 

Method of solving the problem: 

THERMINATOR uses the particle data tables [23] in the universal input form 
used by the SHARE [21] package. The user decides for the thermal parameters 
and the preferred expansion model. The optimum thermal parameters may 
be taken, e.g., as those obtained with the help of SHARE [21] or THERMUS [22]. 
At the moment there are two different expansion models implemented in the 
code: the model of Ref. [16], based on the so-called Buda-Lund [24] param- 
eterization, and the Blast- Wave model [25,26]. The positions and velocities 
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of the particles are randomly generated on the hypersurface according to 
the statistical (Bose-Einstein of Fermi-Dirac) distribution factors. All par- 
ticles, stable and unstable, are included. The particles move along classical 
trajectories from their initial positions, with velocities composed of the ther- 
mal motion and the collective expansion of the system. Stable particles just 
stream freely, while the resonances decay after some (randomly generated) 
time, which is controled by the particle's lifetime. The decays are two-body 
or three-body, and their implementation involves simple kinematic formulas. 
The decays can proceed in cascades, down to the stage where only stable 
particles are present. All particles have tags indicating their parent. The 
secondary rescatterings are not considered in this approach. Full history of 
the event is stored in an output file, allowing for a detailed examination of 
the space-time evolutions and the calculation of the transverse-momentum 
spectra. 

Purpose: 

The ongoing analyses of the SPS and the RHIC data as well as the future 
heavy-ion program at LHC will certainly benefit from THERMINATOR as a 
tool for generating events in a simple statistical model. The Monte Carlo 
code written in C++ and using the standard ROOT [1] environment can be 
easily adapted to purposes directly linked to experimental data analyses. The 
space-time tracking capability will allow, in the framework of the statistical 
approach, to better understand the physics of relativistic heavy-ion collisions. 
THERMINATOR calculates the particle spectra and other observables related to 
the space-time evolution of the system. It provides a C++ framework which 
may be easily developed for detailed analyses of more involved observables 
such as, e.g., correlation functions or HBT radii. 

Computation time survey and storage requirements: 

The generation of 500 events from scratch takes about 1 hour 15 minutes on a 
PC with Athlon-Barthon 2.5 GHz under Red Hat Fedora 3. Each subsequent 
500 events take about 1 hour. To store 500 events about 1.1 GB disk storage 
is needed, depending on the kinematic range. After converting the output to 
the ROOT TTree format, 900 MB may be freed. 

Accessibility: 

http: //hirg. if .pw. edu.pl/en/therminator 
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1 Introduction 



Statistical description of heavy-ion experiments has proven to be very suc- 
cessful in understanding and quantitatively describing the soft physics, i.e., 
the regime with transverse-momentum scales not exceeding, say, 2 GeV. The 
simple basic assumption that the hadronic matter reaches thermal equilib- 
rium and undergoes rapid expansion [27-30] leads to a remarkably good 
description of the ratios of particle abundances measured in heavy-ion ex- 
periments, in particular at the high energies of the SPS {^/s]^=S-17 GeV) 
and RHIC (reaching ^/snn—'^OO GeV). The results for abundances may be 
found in [3-15]. Other details of different statistical approaches may be found 
in [24-26,31-37]. 

Crucial for the success of the statistical approach is the complete inclusion 
of hadronic resonances, whose number grows rapidly according to the Hage- 
dorn hypothesis [30,38-42]. Numerous programs exist to analyze the ratios 
of particle abundances, with SHARE [21] and THERMUS [22] publicly available. 
The abundances are observables that are sensitive to the thermal parameters 
of the fireball created in the collision: temperature, chemical potentials, and 
fugacity factors. The abundances are basically insensitive to the expansion^, 
which means that it only matters how many particles are formed and it is 
irrelevant how fast they move. 

The situation is different for the description of the momentum spectra, 
where the collective motion of the expansion adds up to the thermal one 
generated in the local element of the fireball. Thus, the spectra are sensitive 
not only to the thermal parameters, but also to the geometry of the fire- 
ball and velocity of expansion. Surprisingly, very simple hypotheses for the 
geometry and expansion [16] allow for a good description of the RHIC data 
for the transverse-momentum spectra of identified particles, both in the non- 
strange and strange sectors [43] . The approach of Ref . [16] is based on the 
simplifying assumption that the so-called chemical and thermal freeze-out 
occur simultaneously thus we deal with a single- freeze- out model. Similar 
hypothesis is otherwise known as the sudden hadronization [44]. We follow 
this paradigm in THERMINATOR and throughout this paper ^. 

^Some dependence may come for boost non-invariant systems from the detector accep- 
tance, however in the most studied mid-rapidity region at RHIC and at SPS this is a tiny 
effect. 

^Since the program is an event generator, it would not be very complicated to add an 
"afterburner" that would scatter the particles elastically after they have frozen chemically, 
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THERMINATOR computes the transverse-momentum spectra of identified 
particles in the single-frcczc-out model. The user has an option of choosing 
two variants of the geometry/expansion models: the Buda-Lund type ex- 
pansion [24, 33] and the Blast- Wave type expansion [25, 26] (both models of 
expansion were used with a similar success in the original single-freeze-out 
model [16]). Other expansion models may be looked at with the need of 
only minor modifications to the code. As already mentioned, the key ingre- 
dient in the thermal approach are the hadronic resonances, which must be 
included in a complete way, together with their decay channels and branching 
ratios. THERMINATOR uses the same universal input files with the information 
from the Particle Data Tables [23] as the SHARE package [21] co-authored 
by two of us (WB,WF). The basic functionality of THERMINATOR is the cal- 
culation of the momentum spectra and, correspondingly, the analysis of the 
geometry and expansion. It allows for the optimum determination of the 
transverse flow in the assumed model of expansion. This may be very useful 
for SPS, RHIC, and the future heavy- ion program at LHC. However, the 
program is a versatile Monte Carlo event generator, and as such it may be 
used with slight modifications to analyze other observables, such as HBT 
correlation radii [33,45-48], correlations of non-identical particles, or bal- 
ance functions. With simple modifications amounting to a departure from 
the presently-implemented boost and azimuthal symmetries, one may look 
at the dependence of parameters and observables on rapidity, or at measures 
of the azimuthal asymmetry, such as the coefficient of the elliptic flow V2- 
We want to emphasize that THERMINATOR uses the ROOT environment, thus 
it may be conveniently linked to the standard software used for the data 
analysis. 



2 Basics of the single-freeze-out model 

The basic assumptions of the single-freeze-out model are as follows (for a 
detailed discussion see Refs. [16,37,43]): 

1. At some point of the space-time evolution of the fireball the local ther- 
mal and chemical equilibrium is reached. The local particle phase-space 
densities follow the Fermi-Dirac or Bose-Einstein statistical distribu- 
tions. 

thus differentiating between the two freeze-outs. 
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2. Freeze-out occurs at universal values of the thermodynamic parame- 
ters. We have temperature T and three chemical potentials: baryon, 
fiB, strange, ps, and related to the third component of isospin, /i/3. 
Typically, the potentials ps and /i/3 are set from the condition of the 
overall strangeness neutrality and from the assumption that the ratio 
of the electric to baryon charge density equals the initial ratio, Z/A, 
where Z and A are the atomic and mass numbers of the colliding nuclei. 
The values of the thermodynamic parameters are obtained from fits to 
the measured ratios of the particle abundances. The publicly-available 
SHARE [21] or THERMUS [22] codes may be used here. 

3. The freeze-out occurs on a boost-invariant cylindrically symmetric hy- 
persurface, discussed in detail below. The velocity field of expansion is 
chosen appropriately, providing the longitudinal and transverse flow to 
the system. The particles generated at freeze-out are termed primor- 
dial. 

4. The later evolution consists of decays of resonances which may proceed 
in cascades. All resonances from the Particle Data Tables [23] are 
incorporated. 

The basic physics of the model is depicted in Fig. [T] The popular choices 
of the hypersurface of freeze-out and the velocity field of expansion are dis- 
cussed in detail in Ref. [49]. The boost-invariant formalism has been de- 
veloped by Schnedermann, SoUfrank, and Heinz [25]. The ansatz for the 
boost-invariant, cylindrically symmetric freeze-out hypersurface has the form 



x'' = it,x,y,z) = (r(C)cosh a||,p(C)cos0,p(C)sin0,r(C)sinh ay) . (1) 

The parameter ay is the space-time rapidity. The boost-invariant freeze-out 
hypersurface is completely defined by the freeze-out curve parameterized by 
(. This curve defines the freeze-out times of the cylindrical shells with the 
radius p. The volume element of the freeze-out hypersurface is obtained from 
the equation 

,0^ dx" dx^ dx'^ ,7/1, 

ciS, = .,„,, — — — rfarf/?rf7, (2) 
which with the parameterization yields 

d^Y^^ = I _^cosh an, ^ COS0, — sin0, ^sinh a\\] p(Qf(Qd(^da\\d(j). (3) 
\dC dC dC dC 7 
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Figure 1: The cascade decay in the single- freeze-out model. An unstable 
particle is formed at some position on the hypersurface with the momen- 
tum composed of the thermal motion and the flow of the medium. Then the 
particle travels for some time depending on its lifetime and decays. The prod- 
ucts, if resonances themselves, may decay further, until the stable particles 
are formed. 

Similarly to Eq. (P) the boost-invariant four-velocity field has the structure 

= coshQ;^(^) (coshay, tanha_L(^) COS0, tanha_L((^) sin0, sinh ay) . (4) 

We note that the longitudinal flow is Vz = tanhay = z/t (as in the one- 
dimensional Bjorken model [50]), while the transverse flow at z = has the 
form Vr = tanha_L((^). 

With the standard parameterization of the particle four-momentum in 
terms of rapidity y, azimuthal angle ip, and the transverse mass m^, 

p'^ = {m±coshy, p±_ cos ip, p± sin p, m±smhy) , (5) 

we find 

p ■ u = m^cosh(Q;^)cosh(a|| — y) — p_Lsinh(Q;^) cos(0 — p>), (6) 

and 

p(C)f(C)rfCrfa||rf0. (7) 



£j:-p 



m±cosh.{y - a^)^ - p± cos(0 ~ V^)^ 
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The Cooper-Prye [51] formalism yields the following expression for the 
momentum density of a given species of particles 



dN 



2n 



dyd^PA 



(27r) 



da\\ j dC p(C)f(C) 



X 



m_LCosh(a|| - y)^ - p± cos(0 ~ V^)^ 



X {exp [/5mj^cosh(a^((^))cosh(a|| —y) — 

/3p_Lsinh(Q;_L(C)) cos(0 -if)- P/i] ± 1}"^ , 



(8) 



where (5 — 1/T and /x = -B/xb + Sixs + hl^^h- The expression in curly brack- 
ets is the Bose-Einstein or Fermi-Dirac distribution, with the corresponding 
minus or plus signs. THERMINATOR implements two models of the freeze-out 
hypersurface and velocity field, described in the following subsections. 

2.1 The Cracow single- freeze-out model 

We term the model [16], that has been used very successfully in the descrip- 
tion of the RHIC data, the Cracow single-freeze-out model. It is defined by 
the equations 



f = TCOsh(Q;x), p = Tsinh(Q;x), r = const. 



(9) 



and m'' = x'''/r. This parameterization yields the following distribution of 
the primordial particles 



dN 



1 



dydipp±dp±da\\d(f)pdp {2tt)^ 
X < exp 



m±y^T^~+^ cosh.{a\\ — y) — p_Lpcos(0 — (p) 



o o 
/3m_L \/ 1 + ^cosh(Q;|| - y) - cos(0 -</?)- 



± 1 



(10) 



The variable p is limited from above by the maximum transverse size pn 
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2.2 Blast- Wave model 



The other considered model of the fire-ball expansion is the Blast- Wave 
model ^. In this case one specifies the freeze-out hypersurface and the velocity 
profile of the flow by the conditions: 



T = T = const, Vr = tanha_|_(C) = const. 
Then Eq. (jH)) takes the form 
dN 



T 



dyd(pp±dp±da\\d(f)pdp (27r)^ 



m_|_cosh(a|| — y) 



X < exp 



'1-v: 



cosh(a|l — y) ^ cos(0 — (f) — /3fi 



-1 



± 1 



(12) 



Again, the variable p is limited from above by the maximum transverse size 

Pmax- 



2.3 Decays of resonances 

We incorporate all the four-**** and three-*** resonances. Following the 
scheme implemented in SHARE [21] we exclude all single-* resonances, and 
practically all double-** resonances listed in the Particle Data Tables [23]. 
The reader may use his preferences by modifying the input files. The sta- 
ble particles and resonances are populated according to the formulas (fTUj) or 
(|12p. All particles have their positions and momenta specified. Then, each 
resonance decays after some time, controled by its lifetime l/F. In the rest 
frame of the resonance the decay at time r occurs with the probability density 
r exp(— Fr). The decays are two-body or three-body, and the physical values 
of the branching ratios are taken from the Particle Data Tables. Heavy res- 
onances may decay in cascades. The kinematics of the decay is implemented 
in the rest frame of the decaying particle, and then the system is boosted 
to the reference frame in which the freeze-out hypersurface of Eqs. (jlUj) and 
()12|) is defined. For collisions of symmetric nuclei at RHIC it corresponds to 
the center-of-mass (CMS) frame. 

■^We should stress that the meaning of the blast-wave model here is different from the 
popular usage in fits to the heavy-ion data, since we include the decays of resonances. 
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The particles are distinguished by isospin, hence the branching ratios for 
strong decays incorporate the isospin Clebsch-Gordan factors. The method 
of computing the branching ratios, the pohcy of assigning the branching 
ratios to channels labeled as seen or dominant, and other details concerning 
the decays of particles are incorporated in the same fashion as in SHARE [21] 
where the reader is referred to for more details. 



3 Description of THERMINATOR 

THERMINATOR is written in C++ and uses ROOT - An object-oriented data anal- 
ysis framework [1]. A more sophisticated usage of THERMINATOR at the level 
of analysing the output requires a basic knowledge of ROOT from the user. 

The block structure of THERMINATOR is shown in Fig. |21 Information on 
particle properties as well as the THERMINATOR parameters are read from 
the (modifiable) input files. Two programs, term_events and term_tree 
generate a ROOT tree which can be used within ROOT to generate physical 
results. The blocks of Fig. El are described in detail in the following. 



3.1 Particle data files (particles . data and decays . data) 

The structure of input files containing the information on the particle proper- 
ties and their decays is the same as in SHARE [21], where the reader is referred 
to for a more detailed description. Here we just give the basic information. 

The input file particles . data contains the information on the properties 
of particles such as mass, width, spin, isospin, the quark contents, and the 
Monte Carlo identification number. The file has the following format: 

name mass width spin I 13 q s aq as c ac MC 

where: 



name 


the particle label used in the program 


mass 


mass in GeV, 


width 


width in GeV, 


spin 


spin. 


I 


isospin. 


13 


3rd component of isospin. 


q, s 


number of light / strange quarks. 
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particles. data 
decays. data 
(particle properties) 



therminator.in 
(program parameters) 




therm_events 
(generate events) 



events, out 
(ascii event file) 



fintegrandnnn 
fiTiultiplicity_nnn 
(auxilliary files) 



thermtree 
(convert to TTree) 






ROOT macros 


^ 


generating figures 



physical results, figures, ... 



Figure 2: The block structure of THERMINATOR. 



aq, as number of light / strange antiquarks, 
c, ac number of charm/anticharm quarks, 

MC particle's identification number, where available corresponding to 
the standard Monte Carlo particle identification convention [23]. 
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The file decays. data contains the information on particle decays in the 
format: 

Nameparent Namedaughtori Namedaughtor2 Namedaughters^ BR C-G?(0/1) 

where BR denotes the branching ratio of the decay and C— G refers to 
whether the branching ratio should be multiplied by a Clebsch-Gordan co- 
efficients (0: no, 1: yes). Normally, we use 1 for the two-body decays and 
for the three-body decays, where the full (isospin-dependent) branching 
ratio BR is provided. The entry BR may be used to control the feed-down 
from the weak decays, in particular setting it to would switch off the decay 
channel. 

3.2 Input parameter file (therminator . in) 

The execution of the THERMINATOR program is mainly controled by the param- 
eter input file, whose default name is therminator . in. Different parameter 
filename can be specified by the user at the command line as a parameter 
to the program. The file contains the following information, which must be 
filled by the user: 

NumberOfEvents - the number of events to be generated (default: 500), 
Randomize - if set to 1, the random number generator seed is initialized 
with the help of the current time, producing different seeds for every 
run of the program; if set to the seeds (and therefore all the results) 
will be the same for each run of the program (default: 1), 
InputDirSHARE - location of the SHARE input files (default: ../share). 
Event OutputFile - name of the ASCII file containing the generated events 

(default: event. out), 
FreezeOutModel - selects the version of the freeze-out model used, the 
"Cracow" model (see Sect. 12. 1|) (value: SinglePreezeOut) as well as 
"BlastWave" model (see Sect. 12. 2j) (value: BlastWaveVT) are available 
(default: SinglePreezeOut) , 

BWVt - a parameter specific to the Blast- Wave model - the common ve- 
locity (default: 0.55) {vr from 

Tau, RhoMax - geometric parameters in units of fm (default: 9.74, 7.74) 
(r from ^ or ^ and pmax from Sec. OorQ, 

''^Appears for the three-body decays only. 
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Temperature, Miul, MiuS, MiuB - thermodynamic parameters in units 
of GeV (default: 0.1656, -0.0009, 0.0069, 0.0285) (T from (fTUI) or (fT^ 
and /i/3, fiB, and fis from Sec. EI), 

AlphaRange, RapidityRange - integration ranges for longitudinal vari- 
ables (space-time rapidity and rapidity) (default: 8.0 (-4.0, 4.0), 4.0 
(-2.0,2.0) ) («!! and y from (HHD or dHD), 

NumberOflntegrateSamples - number of samples used in the determina- 
tion of multiplicity and maximum of the integrand (default: 1000000). 

The parameter file is properly commented for the ease of use. 

3.3 Generation of the particle distribution 

The generation of particle distributions proceeds in three main steps, two of 
which are performed only once per given parameter set. 

In order to generate particles through a Monte Carlo method, the max- 
imum value of the distributions on the right-hand-side of Eq. ()10|) or ()12p 
must be known. It is found through a simple numerical procedure. A sample 
of particles is generated and the values of the distributions are calculated 
for each of them. The maximum value obtained is taken as the maximum 
of the requested distribution for the considered particle type. The authors 
have checked, by studying the resulting distributions, that if the size of the 
sample is large enough, the method provides a good estimate of the maxi- 
mum value. The maximum value depends, in principle, on the particle type 
and values of parameters, but does not change from event to event. There- 
fore the value is calculated once for each particle type and stored in the file 
f integrandmax_nnn, where nnn stands for a unique identifier for a given pa- 
rameter set. Subsequent generations of events with the same parameters do 
not require the regeneration of the maximum values; they are read from the 
f integrandmax_nnn file instead. This way the computation time is saved. 

In order to generate events, a multiplicity of each particle type must be 
known. The average multiplicity per event can be calculated in a straight- 
forward manner by numerically integrating the distribution functions in the 
given integration ranges (determined by the model parameters). This pro- 
cedure must also be done only once per parameter set. The average mul- 
tiplicity of each particle is stored in the working directory in a file named 
fmultiplicity_nnn, where nnn is the same unique identifier, as the one 
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used for the maxima discussed above. The multiphcities are then read in by 
subsequent generations, saving the time of calculation. 

In the next step the program proceeds to its ultimate goal of generating 
events, i.e., the data sets containing full information about produced par- 
ticles, the history of their evolution, resonance decays, etc. Each event is 
generated separately. First, the multiplicities of each particle type are gen- 
erated as random numbers from a Poissonian distribution, with the mean 
being the average particle multiplicity determined earlier. Then the program 
proceeds to generate particles, sequentially from the heaviest to the lightest 
particle type. The procedure is, in essence, the generation of the set of six 
random numbers (the magnitude of the transverse momentum pt-, the az- 
imuthal angle of the transverse momentum yj, the rapidity y associated with 
the particle's momentum, the space-time rapidity a\\ associated with the par- 
ticle's location, the transverse distance of the particle p, the azimuthal angle 
of the particle's location 0), distributed according to the formulae (fTUj) or 
(fT^ . depending on the selected source model. The integration ranges are 
determined by the model parameters and the user's input. The event gener- 
ation procedure is a standard von Neumann method of rejection/acceptance 
of the given set of numbers based on the randomly generated test value dis- 
tributed uniformly between and /max, where /max is the maximum value of 
the distribution determined earlier. The accepted sets of numbers are stored 
in memory as representing the actual particles. The procedure ends when 
the determined number of particles of each particle type is generated. At this 
point all the primordial particles, stable and resonances, have been generated 
and stored in the event. 

3.4 Implementation of two-body decays 

The next step in the event generation is the simulation of the decays of 
unstable particles. A particle is treated as unstable if it has non-zero width 
in the particles . data input file. If the user wishes to prevent a decay of 
a particular resonance, he should modify the input file and set the width 
of a given particle to 0. Alternatively, decay channels may be switched-off 
selectively by setting the corresponding branching ratio in decays. data to 
zero. 

The decays proceed sequentially from the heaviest to the lightest particle. 
The decay daughters are immediately added to the set of particles in the 
current event, hence they may decay in the subsequent steps. Most particles 
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have several decay channels. In each decay one of them is selected randomly 
with appropriate probability corresponding to the branching ratio from the 
decays .data input file. It may happen that the selected decay channel is of 
a sub-threshold type (average mass of the parent particle is smaller than the 
sum of masses of the daughter particles), which is not yet implemented in 
THERMINATOR. In such a case, the particle does not decay and the program 
proceeds to the next particle. 

The selected decay channel can be two- or three-body. Both are imple- 
mented in THERMINATOR and treated on equal footing. We turn now to the 
discussion of the two-body decays; the specifics of the three-body decays are 
described in the next subsection. 

First the lifetime Ar of the decaying particle of mass M, mowing with 
the four- momentum p'^, is generated randomly according to the exponential 
decay law, exp(— FAr). Once the lifetime Ar of the particle is known, the 
space-time point of its decay can be calculated, by taking its original space- 
time position a;Qj.igjj^ (determined from the values r, p, ay, and 0), and adding 
to it the space-time distance traveled by the particle, 

^decay ~ ^origin + J^^'^' ^^'^'^ 

This calculation is always done in the CMS reference frame of the colliding 
system in the Cartesian coordinates. Note that Eq. ()13p guarantees that the 
space-time interval between the origin and decay point is equal to Ar. 

The energies Ei, E2, and the three- momenta p^, p2 of the daughter par- 
ticles are initially determined in the rest frame of the parent particle from 
the energy-momentum conservation laws, which yield: 



m 



2 

1,2 



2 = (14) 

. ,_, , _ [(M2-(mi+m2)2)(M2-(mi-m2)2)]V2 
IP1I-IP2I- ^ • (15) 

Here mi, m2 are the masses of the daughter particles. The direction of the 
three-momentum of particle 1 is generated randomly and is evenly distributed 
on the sphere in the parent particle's rest frame. The direction of the three- 
momentum of particle 2 follows from the momentum conservation. Having 
generated the decay momenta of both particles in the resonance rest frame, 
we next boost them to the CMS frame with the parent particle velocity. The 
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decay point of the parent particle (|13|) is taken as the origin of the daughter 
particles. 

Together with each daughter particle the reference to the parent particle is 
stored. This enables the full tracing of the cascade decay tree in the analysis 
of the THERMINATOR output. Also, the parent particle is flagged as "decayed" 
which is later indicated in the output file. 

3.5 Implementation of three-body decays 

Following the common practice, we assume that the transition amplitude in 
three-body decays is constant, i.e., not dependent on the momenta. Similarly 
to the case of two-body decays, the momenta of the emitted particles are first 
determined in the parent particle's rest-frame and only at the end are boosted 
to the CMS frame with the parent velocity. The three-momenta of daughter 
particles considered in the parent particle rest frame lie all in one plane. 
Moreover, the energy distributions of the daughter particles in three-body 
decays with constant matrix element are uniform (see for example Ref. [52]), 
therefore we may generate randomly the energies E2 and of the daughter 
particles 2 and 3 in the energy range allowed by the energy conservation: 
7712 < E2 < M — mi — 7773, 771^ < E^ < M — mi — m2. The energy of the 
daughter particle 1 as well as the magnitudes of all three-momenta (|pi|, IP2I, 
and IpsI) follow from the energy and momentum conservation laws, which in 
addition determine the angle between the vectors p2 and ps. 

El = M — E2 — E^, 

M = ^p^ + 2IP2IIP2I cos 6*23 + pl + ml + ^|p2P + m| + ^jpg^ + m^. 

(16) 

In our approach Eq. (fTH|) is used to calculate cos 623 and whenever the con- 
dition I cos 6^23 1 < 1 is satisfied the generated values of the energies and the 
angle itself are accepted. In this case, we rotate the plane spanned by the 
vectors p2 and pa, together with the whole system, by a randomly chosen 
angle 023 around the axis shown by the vector P3. Next, with the help of 
two additional randomly generated Euler angles the orientation of the vector 
P3 is chosen. Finally, the system is boosted with the velocity of the parent 
particle. Obviously, the generated values of |pi|, |p2|, and |p3| for which no 
solution of Eq. (fTHj) exists are rejected. 
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Similarly to the two-body decays, the space-time origin of daughter par- 
ticles is taken as the decay point of the parent. The daughters are automat- 
ically added to the list of particles in the event. 

3.6 Storage of events 

The particles are stored in the computer memory as collections corresponding 
to one event. Once the whole event is generated (all primordial particles are 
produced and all the unstable particles have sequentially decayed), the event 
is written out to the output file and is erased from the memory. 

Each event in the output file consists of an event header (containing the 
number of particles in the event), and the particles themselves. Each particle 
is represented in a single line. It contains: the particle index for this event 
(useful for cross-referencing the decay parents and daughters). Particle Data 
Group identification number, PID, components of the momentum: p^, Py, 
Pz, energy E, mass, m (in GeV), space-time coordinates of creation: x, y, z, 
t (in 1/GeV), the parent index (-1 if the particle is primordial), finally, the 
decay flag (1 if particle has decayed, if not). An example of the output file 
is shown in Fig. IT^ 

The format of the output file is a plain ASCII text, which is a standard 
used by many event generators. It has the advantage of having the complete 
information and being easily readable by the human or a computer program. 
It is, however, not efficient in terms of disk-space usage and data-access time. 
Therefore a simple ROOT macro therm_tree . C is provided, which converts the 
text output file into a ROOT TTree, which is stored in a ROOT file on disk. 
Events are divided into 500-event batches, each batch in a separate file. 

The information stored in ROOT file differs slightly from the one stored in 
the text file; it does not contain the index. The position coordinates x, y, z, 
and t, are now converted to fm (we use c=l). Instead of the parent index 
and decay fiag the TTree contains two PIDs. The first one, "fatherpid", 
is the Particle Data Group PID of the immediate parent particle (or the 
PID of the particle itself, if it is primordial). The second, "rootpid", is the 
PID of the original primordial particle, that is the particle at the beginning 
of the cascade decay tree. This enables us to distinguish between three 
scenarios: 1) the primordial particle has "pid" = "fatherpid" = "rootpid", 
2) the particle coming from the decay of primordial particle has "rootpid" = 
"fatherpid" and its own, different "pid", and finally 3) the particle coming 
from a cascade decay has its own "pid" , a PID of its immediate parent - the 
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"fatherpid" , and the PID of the original primordial resonance - "rootpid" . 

4 Programming structure of THERMINATOR 

4.1 Technology and general design 

THERMINATOR has been designed as an object-oriented program and has been 
written in the C++ language. It also uses ROOT classes to store and manipulate 
its data structures. 

The classes can be divided into three functional blocks. The first one 
deals with reading input files and storing the information in memory for the 
use by other classes. The second block takes care of storing the information 
of the generated particles and events. Finally, the third block contains the 
main logic of the program, as well as most of the mathematical formulas. 
It deals with the overall program execution control (including parsing the 
parameters), multiplicity calculation, particle generation, the decays, as well 
as storing the results. 

4.2 Data storage block 

The data storage functional block consists of several C++ classes. The top- 
most is the ParticleDB, which stores pointers to information about all Par- 
ticleTypes and allows for easy searches by name or by number. It is filled 
with the data from the input files by the Parser class at the beginning of the 
therm_e vents program. 

The basic object for storing the particle type data is the ParticleType 
class. It contains all the characteristics which are common to all the particles 
of the same type: the charge, baryon number, mass, width, name, etc. It also 
contains the pointer to the complete decay table of this particular particle 
type. There are as many instances of the ParticleType class as there are 
valid particle types read from the input file particles .data. 

The DecayTable class is a collection of pointers to the specific decay 
channels. It deals consistently with the branching ratios. Its job is also to 
select a given decay channel, based on the random number that is passed to 
it. 

The DecayChannel is a simple storage class, which contains the identifiers 
of the daughter particle types as well as the values of the branching ratios 
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for particular decay channels. 

4.3 Calculation results storage block 

This block of classes deals with storing the results generated by the therm_events 
program. It consists of two classes. The general one is the Event class. It 
stores pointers to the particles contained in the event. It is responsible for 
storing them sorted according to mass. It also manages the event generation 
process and holds the generated particle multiplicities for this specific event. 

The Particle class stores information on each generated particle. There 
are as many instances of this class at any given moment, as there are parti- 
cles in a current event. Each particle is defined by its particle type (which is 
a pointer to the corresponding instance of the ParticleType class) and char- 
acteristics that vary from particle to particle, that is the momentum p, the 
creation point x as well as information if it comes from a decay and if so, from 
which particle, and if it has decayed already. This class is also responsible 
for formatting the lines in the output file. 

4.4 Main program block 

This block contains the main acting classes of the therm_events program. 
One of them is the main program, which takes care of the command-line 
parsing, instantiating the ParticleDB classes, initiating reading of parameters 
and inputs, and running the appropriate generation and calculation routines. 

The Integrator class contains the most important part of physics im- 
plementation of THERMINATOR. It includes the procedure which implements 
the Monte Carlo integrator of the emission functions. The integrator has 
three modes of operation: search for the maxima of distributions, straight- 
forward integration (for multiplicity calculation), and generation of random 
numbers according to the probability density described by a given emission 
function. This latter procedure is the particle generation of THERMINATOR. 
All of these modes of operation are properly abstract and are not dependent 
on the particular emission function form. Thanks to this, the emission func- 
tions themselves can be relatively easily modified and new ones can be added 
with little effort. 

One more helpful class included in the main program block is the Par- 
ticleDecayer. Its sole purpose is the implementation of the decay process. 
It assigns the pointer to the father particle. It is responsible for selecting 
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(by the Monte Carlo method) the decay channel, performing the necessary 
kinematic and space-time calculation (again using the Monte Carlo proce- 
dures wherever it is appropriate, e.g., in generating randomly the directed 
decay axis). Finally it returns pointers to the newly created Particle class 
instances containing the daughter particles. Both 2-body and 3-body decays 
are implemented in ParticleDecayer. 

5 Testing against existing calculations 

During the development of THERMINATOR, special care was taken to ensure 
that the obtained results are in perfect agreement with those delivered by the 
MathSHARE package [53]. The most sensitive test available was to check the 
shape of the pr-spectra of many different particle species. Since THERMINATOR 
is based on Monte Carlo methods, the agreement could only be judged within 
the uncertainties of such method, which are decreasing with the increase in 
the number of generated events. 

The agreement at the input level was ensured by simply using exactly 
the same files particle . data and decays. data as input to THERMINATOR 
and MathSHARE. In this way, the same number and types of particles were 
used, as well as the same decay channels and branching ratios. The tests 
were done in two steps. First, the py-distributions of primordial particles 
were compared. The results of this comparison are seen in Fig. |21 For all 
considered particle types, the results of quasi-analytic calculations done in 
MathSHARE and histograms obtained with the help of THERMINATOR agree 
within the error expected from the Monte Carlo procedure. The second 
test was the analysis of the p^-spectra, but this time with full resonance 
contribution included. The results for this test can be seen in Fig. EJ Again, 
the agreement reached with a sample of 1000 events is impressive, proving 
that not only the initialization procedures but also the decay procedures yield 
identical results. 

6 Yielding new results 

This Section contains examples of the physical results obtained with the 
help of THERMINATOR. More results and details focusing on physics will be 
presented in future publications. All the simulations were done with the 
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Figure 3: Comparison of the pr spectra of primordial particles obtained 
from THERMINATOR (circles) and MathSHARE [53] (lines) for the Cracow single- 
freeze-out model. The model parameters corresponding to the best fit to the 
most central PHENIX data for Au + Au @ 200 GeV are taken from Ref. [54]. 
The sample contains 1000 events. The errors on the THERMINATOR results are 
statistical. A perfect agreement of the two independent codes is observed. 

parameters taken from Ref. [54], which were obtained by fitting the STAR 
experiment data with the MathSHARE package. These parameters are included 
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Figure 4: Same as Fig. El for the final (primordial + products of decays) 
distributions. A perfect agreement of the two independent codes is visible. 

as default in the input file therminator . in. The generated sample used here 
contains 500 events. 

In Figs. El and El we show the pT-spectra of various particle species. 
Such spectra were first reported in the original papers on the single-freeze- 
out model [16,43]. The user may generate his own spectra for any desired 
particle by means of simply modifying the attached ROOT scripts f igureS . C 
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or f igure6 . C. 

Figure [7| shows the rapidity distributions of pions, kaons, and protons in 
the midrapidity region. One may observe the boost-invariance imposed on 
the model. With relatively simple modifications the program may be used 
to describe realistic rapidity distribution in the full rapidity range, such as 
for example in Ref. [49]. 

THERMINATOR allows for deeper understanding of the physical results, in 
particular of the role played by the resonances. As an example in Fig. |H1 
we present the anatomy of the vr^ transverse momentum spectra. One can 
see that the contribution from the resonances is dominant, with only one 
quarter of all pions being primordial. The resonance effect is important in 
the whole displayed pt range. One can also see that, as expected from the 
kinematics, the contribution from the three-body decay of the omega me- 
son has a steeper distribution than the contributions from two-body decays. 
Analogous anatomy of the proton px spectra is presented in Fig. IHl 

Further plots show some very interesting and novel observables which 
were not accessible from previous programs implementing the model. These 




Figure 5: The Pt spectra of pions (circles), kaons (squares) and protons (tri- 
angles) generated from THERMINATOR for the Cracow single-freeze-out model 
with default parameters from therminator . in. The sample contains 500 
events. The errors are statistical. The figure is generated by the ROOT script 
figures . C. 
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Figure 6: The px spectra of p mesons (filled circles), Kq mesons (squares), 
mesons (up-triangles), A° barions (down-triangles), and VL~ barions (dia- 
monds) generated from THERMINATOR for the Cracow single-freeze-out model 
with default parameters from therminator . in. The sample contains 500 
events. The errors are statistical. The figure is generated by the ROOT script 
f igure6 . C. 

observables refer to the space-time picture of particle production in relativis- 
tic heavy-ion collisions. For instance, in Fig. we show the distribution 
of the space and time coordinates of the emission points of pions. For the 
primordial particles the emission point is the point on the freeze-out hy- 
persurface, while for particles coming from resonances the emission point 
is simply the point where the decay occurs. The top-left panel shows the 
distribution of the transverse coordinates of the pion emission points, r^. 
The primordial pions are denoted by squares, the pions from the p meson by 
up-triangles, and the pions from the uj meson by down-triangles. The dots 
indicate the full yield (from a// resonances). For r-r in the range from zero to 
Pmax = 7.74 fm (the maximum transverse size of the freeze-out hypersurface) 
the observed distributions are flat, reflecting the uniform distributions of the 
primordial particles. For larger values of tt we can see the tails generated 
by the decays of the resonances. Resonances with longer lifetime generate 
longer tails. This is clearly seen in the top-right panel, where the logarithmic 
vertical scale is used. One can also observe that the contribution from the 
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Figure 7: The rapidity distributions of pions (circles), kaons (squares) and 
protons (triangles) generated from THERMINATOR with the default input file. 
The closed symbols show the rapidity distribution of primordial particles. 
The open symbols show the distributions of the final particles. The sample 
contains 500 events. The boost-invariance assumed in the present model is 
clearly seen. The figure is generated by the ROOT script f igure7 . C. 



long living uj meson saturates the full contribution at tt > 20 fm. 

The two lower panels of Fig. show the analogous analysis with regard 
to the emission time t. The emission starts at t = r = 9.74 fm and peaks 
around t = 14 fm. This reflects the geometry of the freeze-out hypersurface 
assumed in the Cracow model, Eq. ©, as well as the finite lifetime of the 
resonances. 

THERMINATOR may be used to analyze various correlations. As an ex- 
ample. Fig. ^2 shows the space-momentum correlation encountered in rela- 
tivistic heavy-ion collisions. One can see that the x-coordinate of the pion 
emission point is correlated with the a;-coordinate of its momentum. This is 
a characteristic feature of the systems exhibiting the transverse flow. 

The plots shown in this paper are only a few examples of the variety of 
physical results that can be obtained with THERMINATOR. Many novel anal- 
yses may be performed, including a detailed study of the resonance decay 
contribution to particle spectra, a comprehensive study of two-particle cor- 
relations, such as the HBT correlations of pions, kaons, and protons. One 
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Figure 8: The anatomy of resonance contribution to the vr"*" transverse- 
momentum spectra. The plot shows absolute contribution to the overall 
TT"*" spectra (circles) from: primordial pions (squares), pions from the p° 
and decays (up-triangles), and pions from the uj decays (down-triangles). 
The sample contains 500 events. The figure is generated by the ROOT script 
figures. C. 



may investigate also the correlations of non-identical particles, this way get- 
ting the information about relative distances between their emission points. 
All such tasks may be achieved with a complete inclusion of resonance decays. 

An extension of the freeze-out geometry of the model to a cylindrically 
non-symmetric case, such as emerging in non-central collisions, would allow 
for studies of the elliptic flow. As already mentioned, an analogous extension 
to boost non-invariant geometry would prepare ground for studies of the 
dependence of observables on the rapidity variable y. Such extensions are 
planned by the authors in the future. 



7 Installation 

THERMINATOR is distributed in a form of a .tar.gz archive containing the C++ 
source and the data files. Running under ROOT ver. 4. Ox and compiler gcc 
3.x has been tested. In order to run THERMINATOR one needs: 

• C++ compiler 
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Figure 9: The anatomy of resonance contribution to the proton pjn spectra. 
The plot shows absolute contribution to the overall proton spectra (closed 
circles) from: primordial protons (squares), protons from the A° decays (up- 
triangles), protons from the S hyperon decays (down-triangles), and protons 
from the A resonance decays (open circles). The sample contains 500 events. 
The figure is generated by the ROOT script f igure9.C. 

• ROOT libraries and include files 

After unpacking the .tar.gz archive in a separate directory, one compiles 
THERMINATOR by issuing the make command in this directory. The executable 
of THERMINATOR - therni_events is created in the same directory, as well as 
the therni_tree converting program. The SHARE input files are placed in the 
share subdirectory. 

8 Organization of the C++ package 

8.1 Directory tree 
source code (main directory): 

therm_events.cxx - main program 
ParticleDB.cxx, ParticleDB.h - ParticleDB class 
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Figure 10: Space-time picture of the vr"*" emission process. The plots in the 
upper row show spatial distribution of the pion emission points, dN/ {rTdrT). 
The symbol tt denotes the transverse radial coordinate of the points where 
the pions are born. The lower plots show the distribution of the pion emission 
times. Plots on the left are with linear vertical scale, plots on the right are 
copies with the logarithmic scale. The plots show absolute contribution to the 
overall distributions (circles) from the primordial pions (squares), the pions 
from the p° and decays (up-triangles), and the pions from the uj decays 
(down-triangles). The sample contains 500 events. The figure is generated 
by the ROOT script figurel0.eps, see the text for more details. 
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Figure 11: The space-momentum correlation in THERMINATOR. The plot 
shows the correlation between spatial and momentum x-coordinate for the 
pions. The sample contains 500 events. The figure is generated by the ROOT 
script figure 11 . eps 

ParticleType.cxx, ParticleType.h - ParticleType class 

ReadPar.cxx, ReadPar.h - ReadPar class 

Particle. cxx, Particle. h - Particle class 

Decay Table. cxx, Decay Table. h - DecayTable class 

DecayChannel.cxx, DecayChannel.h - DecayChannel class 

Event. cxx, Event. h - Event class 

ParticleDecayer.cxx, ParticleDecayer.h - ParticleDecayer class 
Integrator. cxx, Integrator. h - Integrator class 
Parser. cxx. Parser. h - Parser class 
THGlobal.h - global helper definitions 
Makefile - The makefile 
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therm_tree.C - source of the text-to-TTree conversion macro 
Data files (share subdirectory): 

particles. data - Particle data file 

decays. data - Particle decays data file 
Input files (main directory): 

therminator.in - Default input data file 

9 How to run THERMINATOR 

It is assumed that THERMINATOR has been successfully compiled and the 
therm_events executable is in the user's search path. The user must also 
know the location of the SHARE input files particles . data and decays . data. 

9.1 Instructions for the user 

First, the user must examine/edit the therminator.in file. If the name of 
the file is changed, it then must be supplied as a command-line parameter 
to the therm_events program. The meaning of the various parameters is 
described in section The input file is also conveniently commented and 
contains reasonable default values. A default input file is Fig. ^| Providing 
the input file is in the current directory, one can issue the command: 

therm_events 

This runs the program, which generates the maximum integrand values, the 
average multiplicities, and the events themselves. As a result a text output 
file with the default name event . out is created. Then the user can issue the 
command 

therm_tree 

and the output file is converted into the ROOT TFile event 1 . root containing 
the TTree named particles. If the user changes the name of the output file 
from therm_events to <naine> . out, he must supply this name to therm_tree 
as a command-line parameter. The output file from therm_tree will then be 
called <name>l . root. 
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In case where more that 500 events are requested, the output ROOT files are 
divided into n batches of maximum 500 events, with the names <name>n . root. 

9.2 Sample run 

This sample run assumes that the user is in a directory where the default 
therminator . in file shown in Fig. is located and that the SHARE in- 
put files are in the share subdirectory. After issuing the therm_events and 
therm_tree commands the directory should contain the event . out ascii out- 
put file as well as the event 1 .root output file. The format of the event . out 
file is explained in Fig. ^| 

10 Analysis with ROOT 

The great advantage of the organization of THERMINATOR is the compliance 
to the ROOT standard. Having the event 1. root file, it is possible to browse 
the event tree and generate physical results, even the involved ones, in a very 
simple manner. It is assumed that the user knows at least the basics of ROOT. 
The user can start a ROOT interactive session: 

root event 1. root 

and produce some simple plots from the ROOT command line. For example, a 
plot similar to the squares plot of Fig. Elcan be obtained with the command 
(typed in one line) such as 

particles->Draw("sqrt (px*px+py*py) " , "(1.0/ sqrt (px*px+py*py) ) * 
(pid==211&&f atherpid==21 IMabs (0 . 5*log( (e+pz) / (e-pz) ) ) <0 . 5) " ) 

The contribution from the omega resonance (PDG PID = 223) to the pion 
m-r-spectra can be viewed with a different command, 

particles->Draw("sqrt (mass*mass+px*px+py*py) " , 
"(1.0/ sqrt (mass*mass+px*px+py*py) ) * 

(pid==211&&f atherpid==223&&abs (0 . 5*log( (e+pz) / (e-pz) ) ) <0 . 5) " ) 

The correlation between particle's momentum and position (as seen in Fig.ITT] 
can be obtained with 
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particles->Draw("px:x»hxpx(100, -20, 20, 100,-1.5, 1.5) " , 
"pid==211&&abs (0 . 5*log( (e+pz) / (e-pz) ) ) <0 . 5" ) 

The particle identification numbers, useful in this interactive analysis, can 
be read from the particles .data file. 

More elaborated examples, in particular the macros used to generate the 
figures of this paper (figure5.C, etc.) are included in the THERMINATOR 
distribution. To use data from more than one eventn.root file at the same 
time, use the Add method of the ROOT TChain class several times. An example 
of the method usage in shown in the figure macros. 

11 Conclusion 

We hope that THERMINATOR will become a versatile tool allowing for a simple 
implementation of the thermal approach in analyses of heavy-ion experi- 
ments. It may be useful for theorists exploring their ideas, as well as experi- 
mentalists searching for simple and working models which can be confronted 
to the data. The results of the code can be directly interfaced to ROOT - An 
Object-Oriented Data Analysis Framework. This way, THERMINATOR becomes 
fully compatible with the basic modern software tool of particle physics. 

New possibilities offered by the program were illustrated in Sec. 10 The 
big advantage offered by THERMINATOR is the possibility of exploring the 
space-time evolution of the particle system formed in relativistic heavy-ion 
collisions. This allows for many interesting physics studies which should be 
carried out in the future, such as exploration of various freeze-out surfaces, 
calculation of the HBT correlations with full treatment of resonance decays, 
analysis of the freeze-out criteria, investigation of the role of resonances, cal- 
culation of balance functions in various variables, analysis of correlations of 
non-identical particles, etc. Future extensions of the code should incorporate 
automatic fitting procedures of the model parameters to various observables, 
such as the transverse-momentum spectra or HBT radii. Also, the straight- 
forward extension to cylindrically non-symmetric and boost non-invariant 
geometries is planned, which would allow for further interesting studies. 
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### This is the Therminator input file ### 

# Therminator behaviour is controlled trough this file 

# Please fill in all the neccessary values 

# Number of events to generate 
NumberOf Events = 500 

# Start each event with a new random seed taken from current time 
Randomize = 1 

# The type of input tables 

# Now available : 

# SHARE - input files from SHARE 
TableType = SHARE 

# Specify in which directory the SHARE input f ils are 
InputDirSHARE = ./share 

# Name of the file to store the events in 
EventOutputFile = event . out 

# The model of freeze-out to use 

# Now available : 

# SingleFreezeOut 

# BlastWaveVT 

FreezeOutModel = SingleFreezeOut 

# Values of parameters of the BlastWave model 

# used only when this FreezeOutModel is selected 

# Flow velocity 
BWVt =0.55 

# Parameters common to both models 

# proper time at freeze-out [fm] 
Tau - 9.74 

# maximum transverse radius [fm] 
RhoMax = 7.74 

# Temperature [GeV] 
Temperature = 0.1656 

# Chemical potentials for Isospin, Strangeness, Bar ion [GeV] 
Miul = -0.0009 

MiuS = 0.0069 
MiuB = 0.0285 

# Ranges of integration for z-variables 
AlphaRange = 8.0 

RapidityRange = 4.0 

# Number of samples used in multiplicity and max. integrand determination 
NumberOf IntegrateSamples = 1000000 



gure 12: An example (default) input parameter file therminator . in 
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Therminator tex 






























all_part 
































Therminator 


































. 000 0.000 
































-5 . 606520e-01 


-6 . 202420e-02 


-4 . 821753e+00 


5 


424378e+00 


2 


420000e+00 


-3 . 179843e+01 


-6 . 7733566+00 


-1 .7119616+02 


1 


8113426+02 


_^ 








5 . 972534e-01 


9 . 567281e-01 


-5 . 015253e+00 


5 


681655e+00 


2 


420000e+00 


8 .4838436+00 


1 . 632971e+01 


-9 . 527784e+01 


1 


0890816+02 


-1 








-4 . 55825 9e-0 1 


1 . 6384 07e+00 


-4 . 854391e+00 




650503e+0 




33 9000e+00 


-2 . 54 88 016+00 


2 . 970237e+01 


-9 . 342735e+0 1 




982 666+02 








r T,a 


-8 . 738711e-0 1 


-7 . 882234e-01 


-5 . 668153e+00 




210910e+0 




250000e+00 


-1 . 06 91 63e+01 


-2 . 5 925876+01 


-1 . 248 685e+02 




3719726+02 








d 1 5« 

3^1 Q 


2 . 732464e-01 


-1 . 444380e+00 


-6 . 337567e+00 


5 


8839056+00 


2 


250000e+00 


1 . 1580486+01 


-3 .0779406+01 


-1 . 223908e+02 


1 


3603476+02 


-1 








-5 . 84334 6e-0 1 


1 . 04 57 02e+00 


-5 . 34 6095e+00 




9113516+0 




220000e+00 


-2 . 9507 14e+01 


2 . 158 621e+01 


-7 . 053663e+0 1 




3576316+01 










-9 . 628190e-01 


-4 . 519428e-01 


3 . 557912e+00 


4 


3111616+00 


2 


190000e+00 


-2 .3142066+01 


1 . 485339e+01 


3 . 6687216+01 


6 


7428516+01 


-1 






1 Q 


1 . 922 882e+00 


5 . 978578e-01 


-2 . 4 694 14e-01 




98529 6e+0 




1 90000e+00 


3 . 0317 676+01 


9 . 61 9324e+00 


-1 . 7 62363e+01 




1373936+01 








11 a 


-9 . 148169e-02 


3 . 549711e+00 


9 . 887226e+00 


1 


070271e+01 


2 


045000e+00 


-1 .2429226+01 


2 . 574922e+01 


1 . 0250896+02 


1 


1734596+02 


_^ 





n 


ooo 


7 . 078 858e-01 


-4.4 5 87 94e-01 


-2 . 7 87 678e-01 




2169306+0 




034000e+00 


3 . 8850546+00 


-1 . 3114 976+01 


-8 . 9 93361e+0 




2081196+01 


_-|^ 




11 


?9Q 


-1 . 552147e-01 


-6 . 102028e-01 


-1 . 610166e+00 


2 


6695006+00 


2 


034000e+00 


1 .0011386+01 


-1 .7230436+01 


-3 .3031386+01 


6 


2710946+01 


-1 





10 
1 ^ 


5110 


5 . 620 824e-0 1 


2 . 138 908e+00 


7 . 5254 94e-01 




94 850e+0 




030000e+00 


9 . 8124 996+00 


1 . 14 66 98e+0 1 


-2 . 6275886+00 




17 60566+01 






1 


■^91 S 


-9 . 211756e-01 


-3 . 007009e-01 


2 . 683836e+00 


3 


501837e+00 


2 


030000e+00 


-1 . 353343e+01 


1 .2837126+01 


6 .5102736+01 


8 


384 9776+01 


-1 








-1 . 955413e+00 


1 . 7391 35e+00 


4 . 604367e+00 




6700136+0 




025000e+00 


-1 . 4518386+01 


2 . 6 69975e+01 


6 . 14 38396+01 




4515156+01 






1 ^ 
1 ^ 


SQni 


-1 . 054209e+00 


1 . 208663e+00 


2 . 556273e+00 


3 


6341966+00 


2 


025000e+00 


-1 . 6754546+01 


-1 .3732596+00 


6 .3858066+01 


8 


2499726+01 


-1 






91 O 


5 . 184 602e-0 1 


2 .8787 4 3e-01 


3 . 522405e+00 




100635e+0 




014000e+00 


2 .52667 96+01 


1 . 5 6584 9e+01 


8 . 9224 856+01 




625046+02 






IT 


110 


1 . 181179e+00 


9 . 084622e-02 


3 . 740515e+00 


4 


410339e+00 


2 


014000e+00 


1 .2866916+01 


1 . 545997e+01 


8 .2611336+01 


9 


8354836+01 


-1 








-4 . 442523e-01 


- 5 . 671208e-01 


4 . 370206e+00 




8782486+0 




014000e+00 


-1 . 0832 69e+01 


-1 . 2410516+01 


1 . 27194 86+02 




37 4 5686+02 






1 Q 




5 . 382718e-01 


-1 . 607474e-01 


-6 . 128849e-01 


2 


2020356+00 


2 


OlOOOOe+00 


9 .5379976+00 


-4 .2147726+00 


-1 . 132449e+01 


5 


1782786+01 


-1 





^1 




1 . 449930e+00 


-5 . 741769e-02 


-8 . 508298e+00 


8 


8486856+00 


1 


950000e+00 


3 .1667136+01 


1 . 122480e+01 


-2 .3428006+02 


2 


4178586+02 


-1 






CT^ni 


-1 . 050218e+00 


-9 . 938995e-01 


2 . 305767e+00 


3 


348112e+00 




950000e+00 


-1 . 928524e+01 


-2 .2727936+01 


8 .8560426+01 


1 


0571616+02 


-1 




TO 

^1 




7 . 715302e-02 


1 . 037076e+00 


4. 093344 e+ 00 


4 


6518216+00 


1 


950000e+00 


6 . 9281406+00 


3 . 136835e+01 


1 .5112966+02 


1 


5222406+02 


-1 




^■^ 


~^ 


-2 . 715839e+00 


-7 . 832633e-01 


2 . 378036e+00 


4 


1769416+00 


1 


950000e+00 


-1 .5073886+00 


-9 . 1325376+00 


5 .5284336+01 


7 


4742956+01 


-1 






91 Q 


2 . 414304e+00 


-1 . 524237e+00 


3 . 043907e-01 




4 70924e+0 




950000e+00 


2 . 3021196+01 


-6 . 6080436+00 


4 . 60 91106+00 




5130396+01 








~i 1 1 c 

9511^ 


-9 . 890418e-02 


-8 . 215015e-01 


-1 . 080623e+00 


2 


378002e+00 


1 


950000e+00 


-3 . 14 9344e+01 


-8 .5226506+00 


-8 .2717136+00 


5 


9811106+01 


-1 








5 . 574 685e-01 


5 . 428459e-01 


-8 . 42624 le-01 




25367 96+0 




94 OOOOe+00 


1 .2376966+01 


9 . 430638e+00 


-2 . 014 642e+0 1 




5 610066+01 






o-i 


9-^9t" d 


7 . 017213e-01 


2 . 622612e-02 


-4 . 511965e+00 


4 


9613036+00 


1 


940000e+00 


2 .7106856+01 


-1 .8990496+01 


-1 . 117306e+02 


1 


2658446+02 


-1 




OQ 


9991 


-1 . 6369 95e+0 


8 . 970330e-01 


1 . 24 8185e+00 




9 674 88e+0 




94 OOOOe+00 


-2 . 778232e+01 


2 . 5 90974e+01 


2 . 34 81546+01 




662 6046+01 






TO 


1 1 1 1 C 


-4 . 907981e-02 


6 . 029554e-01 


4 . 021908e+00 


^ 


5018456+0 




930000e+00 


-2 . 07 92 926+01 


1 . 505045e+00 


1 . 3155006+02 


^ 


4216476+02 






■in 


111 IS 


8 . 88350 9e-01 


3 . 8 61457e-01 


5 . 42 1343e-01 


2 


22 64 516+0 




930000e+00 


3 . 4722306+01 


1.8333 50e+00 


-7 . 1 98 091e+0 


g 


0871066+01 










-7 . 480544e-02 


-3 . 928761e-01 


2 . 5308 15e+00 




207783e+0 




930000e+00 


1 . 3027 926+01 


-1 . 600 6076+01 


7 . 222514e+01 




9927016+01 










-6 . 811600e-01 


1 . 125209e+00 


5 . 594173e+00 


5 


062157e+00 


1 


930000e+00 


-1 .8115966+01 


1 .1196266+01 


1 .1014026+02 




2259236+02 


-1 






~i 1 1 C 


3 . 802464e-02 


-2 . 691408e+00 


-9 . 004810e+00 


9 


5946156+00 


1 


930000e+00 


1 . 5335926+01 


-2 .8063306+01 


-1 . 099985e+02 




2489426+02 


-1 




■^^ 

zi 


~9991 ^ 


8 . 292163e-02 


5 . 497412e-01 


2 . 990627e+00 


3 


5971306+00 


1 


920000e+00 


1 .7334476+01 


5 . 755888e+00 


2 .0105246+01 


5 


5522746+01 


-1 






~91 1 1 ^ 


3 . 29 9353e-0 1 


6 . 339179e-01 


4 .20101 8e+00 




673935e+0 




920000e+00 


8 . 9147396+00 


3.0163 60e+01 


1 . 7323596+02 




8287836+02 










-1 . 138183e+00 


-1 . 512093e+00 


2 . 400448e-01 


2 


706641e+00 




920000e+00 


-1 .8573586+01 


-3 .1953306+01 


1 .1800736+01 


6 


2846826+01 


-1 




'i 




5 .24 692 9e-01 


1 . 970 8 OOe-01 


-4 . 1271 02e+00 




5841406+0 




915000e+00 


2 .8662446+01 


1 . 635630e+00 


-1 . 081 820e+02 




2236026+02 










-1 . 522511e-01 


-2 . 420183e-02 


1 . 352760e+00 


2 


3496716+00 


1 


915000e+00 


9 .0190176+00 


-2 .3060616+01 


1 . 9992576+01 


5 


8798916+01 


-1 




39 


13216 


-3 . 865036e-01 


4 . 961109e-01 


-2 . 481840e+00 


3 


197228e+00 


1 


915000e+00 


7 .3514366-01 


-1 .1037776+01 


-5 .4354936+01 


7 


4305496+01 


-1 




40 


-13116 


1 . 386959e-01 


1 . 545150e+00 


6 . 135968e+00 


5 


612416e+00 


1 


915000e+00 


1 .3061306+01 


2 . 065295e+01 


1 .2433326+02 


1 


3601606+02 


-1 




41 


-13116 


-3 . 74797 6e-0 1 


8 . 90207 3e-01 


-6 . 935374e+00 




2594 486+0 




915000e+00 


-9 . 45244 4e+00 


2 . 785301e+01 


-1 . 67 938 6e+02 




7751906+02 






42 


-22222 


-1 . 103887e+00 


3 . 696479e-01 


-3 . 788353e+00 


4 


3994236+00 


1 


910000e+00 


-1 .4403306+01 


1 . 936447e+01 


-8 .5543846+01 


1 


0179286+02 


-1 




43 


-22222 


-1.180230e+00 


-2.232448e+00 


4 . 897472e+00 


5 


8318186+00 


1 


910000e+00 


-8 .2715256+00 


-3.3243426+01 


4.4322656+01 


7 


4716236+01 


-1 


1 


44 


22122 


-1.390629e+00 


5.130754e-01 


-5.086694e+00 


5 


6320206+00 


1 


910000e+00 


-2.8591036+01 


2.6607336+00 


-1 .9688706+02 


2 


0502076+02 


-1 


1 


45 


22122 


1 . 354408e+00 


-3.782093e-01 


2.289813e+00 


3 


2967876+00 


1 


9100006+00 


2.5389576+01 


1.5413556+01 


6.7714316+01 


8 


8948936+01 


-1 


1 


46 


-21212 


-2 . 960950e+00 


-9.416599e-01 


2 . 191330e+00 


4 


2548766+00 


1 


9100006+00 


-3.7822616+01 


2 . 0215536+00 


4 . 9446666+01 


7 


9524116+01 


-1 


1 


47 


1116 


1.928480e-02 


1 .591572e-01 


-7 .169393e-01 


2 


0417476+00 


1 


905000e+00 


1 . 6598026+01 


6.2680386+00 


-8 . 6582686+00 


5 


3237426+01 


-1 


1 


48 


1216 


1.118424e+00 


-1 . 135706e-01 


2 . 403640e+00 


3 


2665406+00 


1 


905000e+00 


2.5234936+01 


7 . 3371046+00 


1 .1307576+02 


1 


2617946+02 


-1 


1 


49 


-212 6 


1 . 386378e-01 


3.051965e-02 


-9.013189e-02 


1 


9124076+00 


1 


9050006+00 


1.0808736+01 


-1.4767156+01 


-5 . 9445476+00 


5 


3053836+01 


-1 


1 


50 


2126 


-7.534297e-01 


-6.843172e-01 


-1.020784e+00 


2 


3889276+00 


1 


9050006+00 


-2 . 6289626+01 


9.0694096+00 


-4.4883926+01 


7 


2335606+01 


-1 


1 


51 


-1216 


-2.133514e-01 


3 .874093e-01 


2 . 438567e-01 


1 


9708116+00 


1 


9050006+00 


7 .2196116+00 


2 . 8068106+01 


3.1102896+00 


5 


7394166+01 


-1 


1 


52 


2226 


2 . 417437e-01 


-7 . 920690e-01 


-5.0057606+00 


5 


4196386+00 


1 


9050006+00 


-2.5890586+01 


-1.1886716+01 


-2 . 1074876+02 


2 


1833716+02 


-1 


1 


53 


2226 


1 . 850948e+00 


6.405938e-01 


-4 .7232896-01 


2 


7728126+00 


1 


9050006+00 


3.0970786+01 


-1.4044726+00 


-3.7734926+01 


6 


9495046+01 


-1 


1 


54 


23124 


-5.010354e-01 


-1 . 086540e-01 


-3.1657716+00 


3 


7225066+00 


1 


8900006+00 


-1.7439956+01 


-2.1305466+01 


-5 . 1746726+01 


7 


5683216+01 


-1 


1 


55 


23124 


-1.354372e+00 


6.543934e-01 


-3.5183236-01 


2 


4409926+00 


1 


8900006+00 


-2.3553966+01 


4 . 4464266+00 


-3 . 9005926+01 


6 


7383206+01 


-1 


1 


56 


1233 


7 . 424609e-01 


-3.788798e-01 


-4 . 6745706+00 


5 


1030086+00 


1 


8693006+00 


1 . 9554846+01 


-1.6762716+01 


-1 .8501926+02 


1 


9419326+02 


-1 





57 


1233 


-2 . 862624e-02 


-2 . 604978e+00 


3.7421196+00 


4 


9279276+00 


1 


8693006+00 


-9.3440576+00 


-2 .1808316+01 


1.0034166+02 


1 


1434966+02 


-1 





58 


337 


7.587479e-01 


-1.497535e+00 


-1.4474926+00 


2 


8897856+00 


1 


8540006+00 


1.5530336+01 


-3.0504306+01 


-1.4680596+01 


6 


1900706+01 


-1 


1 


59 


337 


5 . 963221e-01 


4 .532649e-01 


5.7400896+00 


6 


0784036+00 


1 


8540006+00 


6.9916266+00 


-3.0814876+00 


2.2964866+02 


2 


3503476+02 


-1 


1 


60 


337 


8 . 631210e-01 


9.815413e-02 


8 . 9864486-02 


2 


0493916+00 


1 


8540006+00 


3.5200386+01 


-1 .5051326+01 


-2.5861756+00 


6 


3151986+01 


-1 


1 


61 


13125 


-2 . 659152e-02 


-6.175227e-02 


3.5339356+00 


3 


9802166+00 


1 


8300006+00 


-1.0877096+00 


-1.0814116+00 


9 .5939816+01 


1 


0794116+02 


-1 


1 


62 


13126 


1.111045e+00 


-1 . 174850e+00 


5.2616346+00 


5 


8007226+00 


1 


8300006+00 


7.8703556+00 


-6.1238636+00 


1.3016466+02 


1 


3959496+02 


-1 


1 


63 


13324 


-1.281471e+00 


1.453678e+00 


-1.1159796+00 


2 
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